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Abstract 

The ionization of H2 in intense laser pulses is studied by numerical integration of the time-dependent 
Schrodinger equation for a single-active-electron model including the vibrational motion. The electron 
kinetic-energy spectra in high-order above-thieshold ionization ai^e strongly dependent on the vibrational 
quantum number of the created ion. For certain vibrational states, the electron yield in the mid-plateau 
region is strongly enhanced. The effect is attributed to channel closings, which were previously observed in 
atoms by varying the laser intensity. 
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Above-threshold ionization (ATI) of atoms or molecules by intense laser fields stands 

for the absorption of more photons than needed to overcome the ionization threshold. A simple 
analysis of classical electron trajectories shows that electrons rescattering once from the core after 



the initial ionization step attain large final energies up to lOf/p [IsIjOjO], while direct (unscattered) 
electrons have a maximum energy of 2Up. Here, Up denotes the ponderomotive potential. A 
striking phenomenon arises when high-order ATI is studied with respect to its dependence on 
laser intensity. For the rescattering plateau between 2Up and lOf/p, it was found in experiment 



(a lZD and calculations [|8L fl Q, y H Q, 



1411 . that a slight change in laser intensity can lead 



to order-of-magnitude changes in yield for groups of peaks within the plateau. Explanations were 
found in terms of multiphoton resonances with Rydberg states [18i|9i[l5|] and within the framework 



of quantum paths jl 
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1611 . The spectral enhancements can be related to channel closings that 



occur when the (ponderomotively shifted) lowest ATI peak coincides with an effective threshold 



iili. 



In the present work, we investigate the role of molecular vibration in ATI channel closings, 
using the example of the H2 molecule. Due to the additional degree of freedom, an additional en- 
ergy scale is involved in the dynamics. In the case of atoms, channel closings were introduced by 
scanning through a certain range of laser intensities. We show that for molecules, due to the cou- 
pling between electronic and nuclear motion, intrinsic channel closing effects can be observed by 
comparing the electron- spectra for different vibrational states of the created ion, i.e., applying 
only one laser intensity. 

The minimum energy of a free electron in the presence of a laser field with amplitude Fq and 
frequency cu is equal to the ponderomotive potential (atomic units are used unless stated otherwise) 
Up = Fq /Ato^, which is the quiver energy of the free oscillating electron. Therefore, the laser field 
modifies the ionization threshold. In non-resonant n-photon ionization of an atom, the electron 
carries the final kinetic energy 

Eiiri = nu - Ip - Up, (1) 

with the ionization potential Ip and integer n. The minimum number s of photons needed to free 
a bound electron is therefore defined through s = ceil[(/p + Up)/iu]. By scanning through a range 
of laser intensities and thus varying Up, one can let ATI peaks disappear at the beginning of the 
spectrum. Such a channel closing causes build-up of electron probability near the core Isl] and 
therefore leads to significant enhancements of groups of peaks within the rescattering plateau of 
the ATI spectrum. Note however, that the precise intensities at which this effect occurs deviate 



from the estimate based on the simple formula above lioi fllll. 

Dealing with molecules, due to the coupling of electronic and vibrational motion, energy is also 
transferred to the nuclei of the system, leading to the occupation of vibrationally excited states. 
We expect that for a given vibrational state of the created ion, Eq. ([T]) is changed to 

Eki„ = nu;-/p-t/p-AE^ (2) 

where AE^ = E"" — E'^ is the difference in vibrational energy between the vibrationally excited 
state V in question and the vibrational ground state. Note that Jp denotes here the adiabatic, not 



the vertical ionization potential 12711 . 

Our model of the H2 molecule consists of a single active electron, interacting with two pro- 
tons that are screened by a second (inactive) electron. The electronic motion is restricted to the 
polarization direction of the laser field. We mention that the effect of moving nuclei in ATI of ID 



has been studied previously [|l7|l . but not in the context of channel closings. Note that the 



dynamics of two active e 



ID models of H2 [118 



ectrons coupled to the vibrational motion has been treated earlier within 
1911 . but so far it has not been achieved to calculate ATI spectra within that 
approach. In the present work, the molecular alignment is perpendicular, i.e., the polarization 
direction is perpendicular to the nuclear motion. This choice was made to eliminate the dipole 
coupling between the electronic ground and first excited state of the H2''" ion and hence allow for 



high vibrational excitation of the electronic ground state [|20|]. This leads to the Hamiltonian 



H{z,R,t) =--{ — + 



+ V^iR) + Vi,,iz,R) + E{t)z, 

where the operator E(t)z with E(t) = F{t)Fosm{ujt) describes the interaction of the electron 
with the electric field of a linearly polarized laser pulse in length gauge. The quantity F{t) defines 
the pulse shape and Fq is the maximum field amplitude. The electron coordinate and inter-nuclear 
distance are denoted by z and R; fi^ and fi^ denote the reduced masses of the active electron and 
of the two nuclei, respectively. The vibrational motion in the H2^ ion is incorporated in the model 
by inserting the exact Bom-Oppenheimer ground-state potential of H2''', 

Vn{R) = V^S\r). (4) 

This choice reflects the assumption that the second (inactive) electron stays in the ground state at 
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all times. The active electron interacts with the ion via a soft-core potential 

K.t(^,i?) = -^=i_. (5) 

The idea behind this model is that the core, consisting of the proton charges screened by the 
inactive electron, is treated as one, singly charged object. The values of o-(R) are fitted such that 
the ground-state Born-Oppenheimer potential of the H2 model Hamiltonian matches the exact H2 



Bom-Oppenheimer ground-state potential V^{R) known from the literature [12 ih . A similar fitting 
procedure has been used previously to reproduce the Born-Oppenheimer potential of in a ID 



model [I22I] . Since our model does not allow excitation of the second electron, we have excluded 



the Coulomb explosion channel, which has been extensively studied for example in 12311 . 

The propagation of the time-dependent wave function is based on the split-operator method, 
along with 2D Fourier transformations to apply the kinetic-energy-dependent operators as simple 
multiplications in momentum space. The two-dimensional grid (2-spacing 0.36 a.u., i?-spacing 
0.05 a.u.) extends in i?-direction from 0.2 a.u. to 12.95 a.u., in electronic direction from -276.3 a.u. 
to 276.3 a.u., corresponding to 256 and 1536 grid points, respectively. In the electronic dimension. 



the grid is further extended up to \z\ = 2522.7 a.u. using a splitting technique ll24ll : in this outer 
region the 2D wave function \l/out(^, R, i) is expanded into products states, i.e.. 



where ^j{z,t) and Q{R,t) axe so-called canonical basis states or natural orbitals [|25ll . making 
a single summation over product states possible. They are obtained as the eigenstates of the 
one-particle density matrices of the two coordinates z and R, respectively, for those portions of 
the wave function that are transferred to the outer region. The number of expansion terms is 
chosen to keep at least 99.9% of the total probability. No more than four terms were needed in 
each expansion. Within the outer region, the interaction potential Vint{z, R) is replaced by the R- 
independent potential V^ut(-2) = Vint{z, R = 2). The coupling between z and R is thus removed 
in this area, and ID propagations can be applied separately to the functions ^,j{z, t) and C,j{R, t), 
which allows for large grids. This helps to keep the entire probability on the grid in spite of the 
vast electron excursions in strong pulses. It also leads to high resolution of the kinetic-energy 



spectra. In the range of \z 



10 . . . 161 a.u., the interaction is smoothly varied from Vint to Vc 



out- 



This means that unlike ll24ll . our calculation does not involve a discontinuity of the Hamiltonian at 



the boundary between inner and outer grid. 



The pulse shape F{t) is chosen such that the temporal pulse integral vanishes. The pulses 
have a plateau with constant intensity around the middle of the pulse. The leading and trailing 
edge of the pulses are sin^-shaped ramps. The total pulse length counts five cycles, where 1.5 
cycles are used to ramp the pulse on and off, respectively. The plateau therefore extends over two 
cycles. After the end of the pulse, three empty cycles follow to allow for the continuum electrons 
to significantly escape from the Coulomb potential of the ion before the simulation ends. The 
laser wavelength is 800 nm so that we have Jp = 10.0 u. The system is regarded as ionized for 
\z\ > 30 a.u. The precise choice of this value is not important since the results shown in this work 
involve electrons driven at least 500 a.u. from the ion at the end of the simulation. 

Projecting the ionized part \l/ion (see above) of the final wave function onto the different vibra- 
tional states Xv of H2''', 

M^) = j xtiR) ^ion(^, R, tend) dR, (7) 

leads to one ATI spectrum for each vibrational state considered. The spectra are calculated via 
Fourier transformation to (py^p^), followed by rescaling of the modulus square \(pviPz)\'^ from 
momentum to energy. Only the right-going part of the wave function is used in this work. 

According to Eq. Q, the ATI peaks within these spectra are shifted by AE^. Hence, channel 
closings can be observed by variation of v. Since AE"" can easily exceed the photon energy (for 
the H2'^ ion and A = 800 nm, AE"^ ^ 1.0 u;), the channel closing will always take place, no 
matter where exactly the first ATI peak is located. The different vibrational levels play the role 
of the different laser intensities in the atomic case. As soon as the vibrational energy eats up 
enough energy, we expect that the energy spectrum for the corresponding electrons shows the 
characteristic channel closing features known from atoms. 

In Fig. [H kinetic-energy spectra of ATI electrons belonging to ions that are produced in 
the V = A vibrationally excited state are shown. To enhance readability, except for one example 
only the envelopes are plotted. Since the spectra correspond to a single vibrational state of the 
ion, an atom-like ATI spectrum arises for each laser intensity, and an atom-like channel closing 
can be identified. All spectra contain only data from the z > part of the grid because a slight 
difference in peak positions between z < and z > would lead to blurring of the peaks. The 
envelope top reaches its highest value for an intensity of / = 9.70 x lO^'^ W/cm^ (Up = 3.73 to) 
at around 42 eV. Trying to estimate this intensity from in Eq. the term AE"" is kept constant 
while Up is scanned through. Using Eq. we find the channel closing (-Ekin = 0) to be located at 
Up = (k + OAO) u, with integer, non-negative k, where k + 11 = n, corresponding to a higher laser 
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FIG. 1: Envelopes of kinetic-energy spectra of (right-going) ATI electrons leaving the ion in the v = A 
vibrationally excited state. The various laser intensities correspond to ponderomotive potentials between 
Up = 3.13 w and Up = 3.93 uj. 

intensity of approximately I = 1.14 x lO^'* W/cm^ for n = 15 or a lower one of approximately 
/ = 8.82 X 10^^ W/cm^ for n = 14. Note that these values refer to our model Hamiltonian using 
fitted potentials and ignoring small effects such as mass polarization. 

If the laser intensity is kept fixed, but the spectra for several vibrational states of the created 
ion are plotted (see Fig. [2]), an intrinsic channel closing (ICC) appears. In this case, in Eq. 
Up = 3.73 u; is kept fixed and AE'" is scanned from f = to i; = 7. This corresponds to values 
between AE^ = 0.00 uj and AE"^ 1.02 cu. The amount of energy the electron loses in each case 
as a consequence of energy conservation shifts the spectra with respect to each other. The ICC 
is masked by the fact that the different vibrational states are not populated equally (see Fig. [3]). 
Therefore, normalized spectra are plotted in Fig. Oa), where each spectrum has been divided by 
the total yield of the corresponding vibrational state. The v = A spectrum shows the highest yield 
in the middle hump compared to all other spectra within the plot. We attribute this behavior to 
an ICC. The unnormalized spectra in Fig.|2b) show the highest yield of the middle hump already 
at V = 2 due to the suppression of higher vibrational states. The application of Eq. ^ leads to 
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FIG. 2: Envelopes of kinetic-energy spectra of (right-going) ATI electrons produced by a laser pulse of the 
intensity / = 9.696 x 10^'^ W/cm^ (Up = 3.73 uj), plotted for vibrational states of the H2''~ ion from 
V = 0, . . . ,7. In (a), the spectra are divided by the total yield of the corresponding vibrational state (cf. 
Fig. [3]). In (b) the spectra show the correct weighting with respect to their yield. 
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vibrational level 

FIG. 3: Occupation of vibrational states after a laser pulse with an intensity of / = 9.696 x 

(Up = 3.73 iu), split up into contributions from the left and right parts of the grid. The Franck-Condon 

overlap is plotted for comparison. The distribution has been normalized to total probability one . 

AE'" = {k + 0.31) Lo, where k + lA = n, which corresponds to an energy difference close to the 
vibrational state v = 2 (using n = 14). Again, the observed position of the channel closing is 
shifted with respect to the one expected from Eq. Q. 

It should be stressed that the energy difference between two vibrational states is larger than 
0.1 uj for all vibrational states considered in this work. Hence, in contrast to the intensity scanning, 
the transition over a channel closing is not sampled continuously. Similar calculations for D2 are 
work in progress and provide a finer graining due to the closer-lying vibrational states of the 
ion. 

We show in Fig. HI that for a suitable electron-energy window the ICC structure appears in 
the electron yield plotted as a function of vibrational quantum numbers. The distributions are 
shown using either the normalized electron yield from Fig. 12 a), see Fig.Ua), or the unnormalized 
electron yield from Fig. [2tb), see Fig. Ilfb). We use the energy window between 40 and 50 eV, 
corresponding to the dashed lines in Fig. [IJ Clearly the ICC shows up in Fig. Ufa). Although the 
ICC feature is not as evident in the unnormalized distributions of Fig.Hfb), it is clearly visible that 
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FIG. 4: Electronic yield within the energy window 40-50 eV (cf. Fig. [T]), plotted versus vibrational states 
of the H2'^ ion for different laser intensities. In (a), the electron spectra used were divided by the total yield 
of the corresponding vibrational state (cf. Fig. |3]l as in Fig. Oa). In (b), correctly weighted spectra with 
respect to their yield were taken as in Fig.|2l;b). 
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electrons and ions are highly correlated, since the distribution over vibrational states is, within the 
chosen electron-energy window, very different from the general distribution of vibrational states 
for all electrons. See, e.g., the curve for Up = 3.73 lu as compared to Fig.[3l 

To summarize, we found clear signatures of spectral enhancements due to channel closings oc- 
curring in ATI of H2 by scanning through the vibrational states of the created ion. The expla- 
nation of this effect seems straightforward, applying energy conservation to the photon-absorbing 
molecule. Similar to atoms, we find that the effect occurs at intensities/vibrational states slightly 
different from the simple estimate based on the unperturbed ionization potential. We conclude with 
a note on the experimental perspectives. The populations of H2"*' vibrational states after strong- 
field ionization of H2 has been measured in [|20n . but a measurement in coincidence with electrons 
will be difficult. On the other hand, coincidence measurements similar to recent pump-probe ex- 
periments [|26n appear feasible. The goal would be to measure the electron from an ionizing pump 
pulse, together with fragments from probe-pulse-induced Coulomb explosion of H2^. 

This work has been supported by the Deutsche Forschungsgemeinschaft. 



[1] P. Agostini, F. Fabre, G. Mainfray, G. Petite, and N. K. Rahman, Phys. Rev. Lett. 42, 1127 (1979). 

[2] J. H. Eberly, J. Javanainen, and K. Rz^zewski, Phys. Rep. 204, 331 (1991). 

[3] G. G. Paulus, W. Nicklich, H. Xu, P Lambropoulos, and H. Walther, Phys. Rev. Lett. 72, 2851 (1994). 

[4] G. G. Paulus, W. Becker, W. Nicklich, and H. Walther, J. Phys. B 27, L703 (1994). 

[5] G. G. Paulus, W. Becker, and H. Walther, Phys. Rev. A 52, 4043 (1995). 

[6] P Hansch, M. A. Walker, and L. D. Van Woerkom, Phys. Rev. A 55, R2535 (1997). 

[7] M. P. Herdein, P H. Bucksbaum, and H. G. Muller, J. Phys. B 30, L197 (1997). 

[8] H. G. Muller and F C. Kooiman, Phys. Rev. Lett. 81, 1207 (1998). 

[9] H. G. Muller, Phys. Rev. A 60, 1341 (1999). 

[10] G. G. Paulus, F Grasbon, H. Walther, R. Kopold, and W. Becker, Phys. Rev. A 64, 021401(R) (2001). 

[11] R. Kopold, W. Becker, M. Kleber, and G. G. Paulus, J. Phys. B 35, 217 (2002). 

[12] S. V. Popruzhenko, P A. Komeev, S. P Goreslavski, and W. Becker, Phys. Rev. Lett. 89, 23001 (2002). 

[13] J. Wassaf, V. Veniard, R. Taieb, and A. Maquet, Phys. Rev. Lett. 90, 013003 (2003). 

[14] K. Krajewska, 1. 1. Fabrikant, and A. F Stai'ace, Phys. Rev. A 74, 053407 (2006). 

[15] R. M. Potvliege and S. Vucic, Phys. Rev. A 74, 023412 (2006). 



10 



[16] R. Kopold, W. Becker, and M. Kleber, Opt. Commun. 179, 39 (2000). 

[17] A. D. Bandrauk, S. Chelkowski, and I. Kawata, Phys. Rev. A 67, 013407 (2003). 

[18] T. Kieibich, M. Lein, V. Engel, and E. K. U. Gross, Phys. Rev. Lett. 87, 103901 (2001). 

[19] S. Saugout, C. Cornaggia, A. Suzor-Weiner, and E. Charron, Phys. Rev. Lett. 98, 253003 (2007). 

[20] X. Urbain, B. Fabre, E. M. Staicu-Casagrande, N. de Ruette, V. M. Andrianarijaona, J. Jureta, J. H. 

Posthumus, A. Saenz, E. Baldit, and C. Cornaggia, Phys. Rev. Lett. 92, 163004 (2004). 
[21] W. Kolos, K. Szalewicz, and H. J. Monkhorst, J. Chem. Phys. 84, 3278 (1985). 
[22] B. Feuerstein and U. Thumm, Phys. Rev. A 67, 043405 (2003). 

[23] S. Chelkowski, A. D. Bandrauk, A. Staudte, and P B. Corkum, Phys. Rev. A 76, 013405 (2007). 
[24] R. Grobe, S. L. Haan, and J. H. Eberly, Comp. Phys. Comm. 117, 200 (1999). 
[25] P-0. Lowdin, Phys. Rev. 97, 1474 (1955). 

[26] T. Ergler, A. Rudenko, B. Feuerstein, K. Zrost, C. D. Schroter, R. Moshammer, and J. Ullrich, Phys. 

Rev. Lett. 97, 193001 (2006). 
[27] The adiabatic ionization potential is the difference between the total ground state energies of H2~'^ and 

H2 (including nuclear motion). 



11 



